Lines of Fisher's zeros as separatrices for complex renormalization group flows 



O 
(N 



in 



> 

oo 

rn 
O 



Yuzhi Liu and Y. Meurice 
Department of Physics and Astronomy, The University of Iowa, Iowa City, Iowa 52242, USA 

(Dated: January 18, 2013) 

We extend the renormalization group transformation based on the two-lattice matching to the 
complex inverse temperature plane for Dyson's hierarchical Ising model. We consider values of 
the dimensional parameter above, below and exactly at the critical value where the ordered low 
temperature phase becomes impossible for a real positive temperature. We show numerically that, 
as the volume increases, the Fisher's zeros appear to accumulate along lines that separate the flows 
ending on different fixed points. We justify these findings in terms of finite size scaling. We argue 
that the location of the Fisher's zeros at large volume determine the phase diagram in the complex 
plane. We discuss the implications for nontrivial infrared fixed points in lattice gauge theory. 
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Classical field equations suggest the existence of mass- 
less Nambu-Goldstone modes in low temperature spin 
models with continuous symmetries, or massless gauge 
bosons in weakly coupled models with local symmetries. 
However, in some cases, quasiparticle excitations or in- 
stantons destroy the long range-order and generate a 
mass gap j| . This theoretical framework provides a justi- 
fication for the absence of long-range order for the D = 1 
Ising model or the D — 2 nonlinear 0(3) sigma model 
with nearest neighbor interactions, and the D = 3 U(l) 
or D = 4 SU{2) lattice gauge theories with a compact 
Wilson's action. 

For gauge theories, the absence of long-range order 
is associated with confinement. The possibility of al- 
most losing confinement by adding a suitable number of 
fermions, in order to have a "walking" coupling constant 
has been considered for models of electro-weak symme- 
try breaking beyond the standard model This has 
motivated an intense activity in the lattice gauge theory 
community [3, 0] ■ From this point of view, it is crucial 
to decide unambiguously if an unusual infrared behav- 
ior will be present for a given number of fermion fields. 
More generally, proving the (non) existence of a mass gap 
is often important and difficult for condensed matter and 
particle physics models. 

In the renormalization group (RG) approach, confine- 
ment means that the RG flows go uninterrupted from 
the weak coupling region to the strong coupling region 
[3, II] . Recently, the loss of long-range order or the ap- 
pearance of confinement when a parameter is varied has 
been explained in terms of RG fixed points disappear- 
ing in the complex coupling or temperature plane [6|-l9j. 
In this Letter, we provide a direct illustration of this 
mechanism with numerical calculations of the RG fiows 
in the complex inverse temperature (0) plane using the 
two-lattice matchi ng p rocedure 13) U | for Dyson's hier- 
archical model with an Ising measure. The choice 
of this model is justified by the feasibility of difficult nu- 
merical calculations explained below. The complex RG 
flows for three values of D are illustrated in Fig. [T] 



From a practical point of view, our main point is that 
as the volume increases, the zeros of the partition func- 
tion in the /3 plane (Fisher's zeros) accumulate along lines 
that separate the flows ending on different fixed points 
and thus determine the complex phase diagram. 



Complex RG flows for various dimensions 




FIG. 1. RG flows in the complex /? plane for D = 3 (top), 
2 (middle) and 1.7 (bottom) and Fisher's zeros for rimax ~ 3 
and 4. A darker background indicates an ambiguous solution. 
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The Hamiltonian of Dyson's hierarchical model with 
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where -B*^"-' denotes hierarchically nested blocks of size 
2". For details, we refer to a recent review article [l^ . 
The parameter c controls the decay of the interactions 
with the size of the blocks. From the scaling of a free 
massless Gaussian field under a change of the lattice spac- 
ing by a scaling factor b, we can include a dimension D 
through the relation c/4 = b^^~^ . Since the number of 
sites is divided by 2 at every blockspin transformation, we 
have = 2. For D = 3 (c = 2^/'^), and more generally 
for D > 2, the Ising hierarchical model has a second order 
phase transition and has many common features with the 
D = 3 nearest neighbor interaction Ising model, however 
the critical exponent f] is zero. For D — 2 {c — 1), and 
more generally for D < 2 the model has no phase tran- 
sition at finite temperature unlike its Z3 = 2 nearest 
neighbor interaction counterpart. Our numerical study 
will focus on D = 3, 2 and 1.7. The case D = 2 is at the 
boundary and is quite interesting, as are other models 
playing this role (for instance, the D = 2 nonlinear 0(2) 
sigma model with nearest neighbor interactions). 

The partition function is obtained by integrating e~^^ 
with a local measure W^o(0)- It can be calculated iter- 
atively because the block variables are not mixed in H. 
If we call Wn((/)) the unnormalized, unrescaled, probabil- 
ity distribution of (/> in a block of size 2", we obtain the 
(implicitly /3-dependent) recursion relation : 



1 / c \Ji + l J 2 



(2) 



In the following, we consider an Ising measure Wo{4>) = 
5{(tP' — 1) and Eq. © reduces to finite sums. This allows 
an exact calculation of Wn{(f) as a function of (3 as long 
as we have enough computer memory to keep track of 
the terms. At finite volume, the partition function is an 
entire function of the form: 



k 



(3) 



The zeros of Z can be obtained from the intersections of 
the zero level curves for the real and imaginary parts or 
from logarithmic residue methods used in Ref. |15| . 

In Fig. [21 we compare the lowest (closest to the real 
axis) Fisher's zeros for different volumes. The distribu- 
tion of the Fisher's zeros is symmetric about the real axis 
and we only show the zeros in the upper half plane. For 
the D = 2> case, the system has a second order phase 
transition at /3c — 1.179 [l3| and we can check the con- 
sistency of our results with finite size scaling 1^ 17 1. In 
Fig. m the lowest Fisher's zeros accumulate toward (3c 
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FIG. 2. The lowest Fisher's zeros for _D = 3 (filled circle 
ending on /3c), 2 (crosses) and 1.7 (filled square) with rimax 
from 2 to 11. As n^ax increases, the imaginary part decreases 
for D =3 and 2 (top panel). Real and imaginary parts of 
Fisher's zeros for _D = 3 and rimax going from 2 to 11 (bottom 
panel). 



as the volume increases. The departure from a linear be- 
havior is significant and requires subleading corrections. 
The singular part of the free energy can only depend on 
the linear size of the system L through the combination 
KiL^', where Ki are the nonlinear scaling variables. Un- 
der a RG transformation, L L/b and Ki — > U^^Ki. If 
we only keep the relevant variable (yi = l/i') and the 
first irrelevant variable (y2 = — w), the requirement that 
/3i(L) is the (lowest) zero imposes a relation of the form 



KL^/" = A + BL- 

Using the approximation K ~ 
zero, we obtain 



+ OiL-^'^) (4) 
/3i(L) — f3c for the lowest 



Re{(3i{L)) ~-(3c = Re{A)L-^/'' + Re{B)L-^/''-'^ (5) 
Im{Pi{L)) = Im{A)L-^/'' + Im{B)L-^/''-'^ 

Using the known values for /3c, v, and ui [14i | and fitting 
for A and B, we obtain results displayed in Fig. [2] The 
fits are in very good agreement with the data except for 
the smallest volume. 
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Complex RG flows have been studied and discussed for 
a variety of statistical mechanics, many-body and lattice 
models Here we construct the RG flows 

in the complex /3 plane b y using a two-lattice matching 
method inspired by Refs. 13, hj. At level n in the block 
hierarchy, we split the system into two blocks Bi and B2 



each with 2" sites. The observable we considered is 



eB2 



Si <l>y )) 



(6) 



0,n 



This ratio measures how the blocks are correlated. It 
does not require any division by the partition func- 
tion. The field rescaling which needs to be calculated to 
write down a full RG transformation cancels out. Using 
Wn-iifj)) for the two block variables, we obtain 

i?(/3,n)= (7) 
/ d0id(/.2e4(S)"(^i+-^^)Vi02W„-i(./.i)W^„-i(./.2) 

Unlike RG transformations based on decimation for one- 
dimensional Ising models 3, 1^ , the measures obtained 
by iterating Eq. ([2]) are not Ising measures. Instead, we 
rely on the assumption that after enough RG transforma- 
tions, the flows become approximately one-dimensional 
and that we can use /? as a coordinate for the unstable 
direction. The matching condition reads: 



i?(/3,n) = R{P',n~l) 



(8) 



with n as large as possible. Given an initial /?, the nu- 
merical calculation of (3' can be done with the methods 
used to find Fisher's zeros. 

This procedure yields very simple flows for real values 
of /3. They can be described in terms of a discrete Callan- 
Symanzik /3 function defined as: 



Am = 



(9) 



The existence of zeros of A/3 implies the existence of 
fixed points. Figure [3] shows A/3 for I? = 3, 2 and 1.7. It 
is not surprising that there is only one zero at the origin 
for _D = 2 and 1.7 since it has been proven fl3| that there 
is no phase transition at finite temperature. For D — 3, 
another nontrivial zero appears near /3c. As the volume 
increases, the nontrivial zeros approach /3c. 

For D = 2, a low temperature expansion can be used 
to show that at leading order 



A/3(/3) 



(10) 



The lowest order contributions to i?(/3, n) for Z) = 2 
comes from flipping all the spins in one of the blocks 
of size 2"~^. The energy cost is only -1/2 and is n- 
independent. The leading order contribution to A/3 
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FIG. 3. Discrete 13 function for D= 3, 2 and 1.7 



comes from the flrst rt-dependent energy cost. Numeri- 
cally, B3 = 35/32,^4 = 155/128 and B5 = 651/512. The 
general formula is B„ = ((ll-hl/4"-i)/6)-l/2"-i-l/2 
and i?oo =4/3. The exponential decay is in contrast with 
A/3 for D = 2 nearest neighbor 0{N) models in the large 
N limit which has been calculated for the two-lattice 
matching [lO| and a simple rescaling of cutoff in the gap 
equation [9|, where at infinite volume, A/3 approaches 
ln6/(27r) asymptotically. It is possible to find simple 
models for the continuous Callan-Symanzik /3-function 
with continuous parameters that interpolate qualitatively 
among the various behaviors of A/3 in Fig. [3] and that 
have zeros moving in the complex plane when D — 2. 
Functional conjugation methods provide a more sys- 
tematic approach of the relationship. 

In the complex case, it appears in general the match- 
ing equation has more than one solution. For D = 2, 
there is a fixed number of solutions, because the match- 
ing equation reduces to a polynomial equation (the Ek 
in Eq. ([3]) are all rational numbers). In order to resolve 
this ambiguity, we picked the /3' that minimizes |/3 — /3'|. 
Under some circumstances, there is only one /3' close to 
/3 and this seems quite natural but this is not always the 
case. In order to quantify this ambiguity, we defined: 



/(/3) = log 1/3-/3^1 -log 1/3-/311 



(11) 



where f3[ is the closest solution and the second closest. 
In Fig. [H we made a contour plot of f{(3). The darker 
the color, the more ambiguous is the selection of /3'. In 
dark regions, the flow is rather erratic. This reflects the 
existence of several competing solutions rather than some 
intrinsic "chaos" as for the D — 1 Ising model On 
the other hand, the mostly unambiguous flows of Fig. [1] 
show simple patterns. For D = 3, the zeros appear near 
the boundary of the basins of attraction the two stable 
fixed points at and 00. As D decreases, the fixed point 
moves to larger /3 and becomes infinite at D = 2. For 
lower D, a pair of complex conjugated zeros appears. 
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D=2, zeros up to n=ll 
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FIG. 4. Detail of Fig. [T] for D = 2 and Fisher's zeros for 
larger number of sites described in the text. 

Figure [T] indicates that Fisher's zeros for the two vol- 
umes matched appear in regions where the flows are am- 
biguous and at the end of flow lines that neither go to zero 
or other fixed points. The zeros for larger volumes move 
approximately "backward" along the set of flow lines that 
separate the flows going to different fixed points. This 
is illustrated in Fig. 2] for D ~ 2, where all the zeros 
corresponding to volumes up to 2^ fitting in the frame 
and the lowest zeros for volumes up to 2^^ have been 
displayed. This behavior is in approximate agreement 
with the argument 18] that under a RG transformation, 
Fisher's zeros for 2" sites should map into Fisher's ze- 
ros for 2"~^ sites. In view of this argument, the general 
shape of the flows and the absence of transition on the 
real axis, it is plausible that as the volume increases, the 
line extends to infinity. Lines of Fisher's zeros are quite 
common T?] . Surfaces were also observed in the complex 
variable sinh(2/3) [211. 

In summary, we have calculated the complex RG flows 
for three values of D. For D = 3, the Fisher's zeros 
pinch the real axis at /3c and separate the flows going to 
(symmetric phase) and oo (broken symmetry phase). 
For D < 2, it is plausible that the line of Fisher's zero 
goes to infinity with no contact with the real axis. This 
hypothetically infinite line separates the flows going to 
from those curling back to either a complex fixed point or 
infinity. This suggests that the basins of attraction of the 
various fixed points in the complex plane are separated 
by lines of Fisher's zeros. From a practical point of view, 
it is easier to calculate Fisher's zeros than to construct 
complex RG flows. A direct interpretation of complex 
flow would be desirable. In problems involving oscilla- 



tions, the addition of an imaginary part to the energy 
or frequency corresponds to dissipation, skin depth or a 
finite lifetime. A simple example is the electric permit- 
tivity e becoming complex as a result of damping. From 
this point of view, a complex energy density (e/2)i?^ is 
an effective description of the interactions with matter. 
Complex eigenvalues for the transfer matrix were also 
found ^] and could suggest a more complete picture of 
complex RG flows. 
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